Can experiments select the configurational component of excess entropy? 



Silvia Corezzi, Lucia Comez, and Daniele Fioretto 
1 Dipartimento di Fisica and INFM, Universitd di Perugia, 1-06100, Perugia, Italy. 

(Dated: November 15, 2002) 

We introduce an experimental method of assessing the vibrational and configurational compo- 
nents of the excess entropy of a liquid over crystal, based on a joined investigation of dynamic and 
thermodynamic properties as a function of temperature and pressure. We analyze light scattering, 
dielectric, calorimetric and dilatometric measurements of three prototype glass formers, orthoter- 
phenyl, salol, and glycerol. In all cases we find that about 70% of the excess entropy is configurational 
in nature. 

PACS numbers: 



Liquids become glasses as a consequence of the sup- 
pression of long range molecular mobility. The most ev- 
ident feature that accompanies glass formation is a dra- 
matic increase of the structural relaxation time, r a . A 
key for understanding this phenomenon is the Adam- 
Gibbs (AG) theory Q), which states a relationship be- 
tween r a and the configurational entropy S c , 

Here, S c = S melt -S vibr where S melt is the entropy of the 
metastable liquid phase, and S' mbr is the entropy of an 
ideal amorphous-solid phase (ideal glass) in which only 
vibrations are active; C is nearly constant, and To is the 
relaxation time in the high temperature limit. 

The difficult point in testing the predictions — and the 
validity itself — of Eq. ([!]) is the evaluation of the config- 
urational entropy from experimental data. In the past, 
this problem has been overcome by approximating the vi- 
brational heat capacity at low temperatures by the heat 
capacity of an experimentally accessible solid phase, i.e., 
the crystal or the glass Therefore, configurational en- 
tropy is replaced by excess entropy. Notwithstanding, the 
AG equation succeeds in describing experimental results. 
This fact looks like a paradox, since the same evidence 
is also found in computer simulation studies, where the 
actual configurational entropy is calculated. An approx- 
imate proportionality between configurational and total 
excess entropies could be a possible explanation of the 
paradox. From Eq. (|l|) it is evident that relaxation time 
measurements performed as a function of temperature 
alone, cannot help to distinguish between S c and a quan- 
tity proportional to it, as the proportionality constant 
simply renormalizes the value of C. 

Numerical simulations of supercooled water, and ana- 
lytical models for thermodynamics of defect crystals (see 
Ref. ||), seem to support this idea. On the other hand, 
to our knowledge, the only experimental work assessing 
the two separate contributions to excess entropy has been 
made on a very simple material, elemental selenium [Q . 
In this case, the vibrational entropy over crystal is calcu- 
lated under harmonic approximation from the analysis of 
neutron-scattering data, and it is subtracted from S exc to 
determine S c . This result gave rise to controversial con- 



clusions about the method of calculation of S c , and the 
existence of proportionality between S c and S exc (|| 

Here, we enter the debate trying to answer the ques- 
tion: "Can experiments select the configurational com- 
ponent of excess entropy in glass forming systems?" . We 
show that calorimetric, dilatometric, and relaxation mea- 
surements as a function of both temperature T and pres- 
sure P, allow us to extract the configurational fraction of 
the excess entropy. Our method is tested on three differ- 
ent prototype glass formers, orthoterphenyl, salol, and 
glycerol and for relaxation data, obtained by dielectric 
spectroscopy and light scattering. 

The configurational entropy of a system at a given T 
and P can be separated into an isobaric contribution at 
zero pressure, and an isothermal contribution at temper- 
ature T @: 

S C (T, P) = S i c sob (T, 0) + S l c soth (T, P) = 

(2) 

where ACp — C™ elt - C P lbr is the configurational heat 
capacity at zero pressure, A (dV/dT) p = (dV/dT)™ elt — 

(dV/dT) V p hr is the configurational thermal expansivity 
at T, and Tjc is the temperature where the configura- 
tional entropy of a liquid vanishes at P = 0. Substi- 
tuting Eq. (0) in Eq. (fil) gives the pressure extended 
AG equation (PEAG) g7 For evaluation purposes, also 
in this approach one approximates the isobaric part of 
the configurational entropy with the excess entropy, i.e., 
Si sob (T,0) ra S exc (T,Q). Within this approximation, the 
PEAG model is able to describe experimental data fairly 
well in the T-P space |6|, fjj. The cost to be payed is a 
value of (dV/dT)p which m some cases may consid- 
erably depart from the expected value. This anomaly 
is a consequence of the poor approximation S c ra S exc . 
In contrast, in this work we improve the PEAG model, 
we obtain a reasonable value for (dV/dT) V p br , and, more 
interestingly, we determine a realistic estimate of the con- 
figurational entropy. 

The model. — Guided by the idea that an approximate 
proportionality may exist between configurational and 
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FIG. 1: a) Structural relaxation time of OTP from dielectric measurements at different pressures. Open symbols are data taken 
at different pressures Q (P=0.1 MPa (o), 19.6 MPa (<), 39.2 MPa (A), 58.8 MPa (0), 78.5 MPa (□)). Crosses are data at 
ambient pressure 0|. d) Structural relaxation time (r) of OTP from photon-correlation measurements at different pressures. 
Data taken from Fytas et al.fuj (P=0.1 MPa (<), 25 MPa (□), 50 MPa (o), 75 MPa (>), 100 MPa (*) 125 MPa (0)). c) 
Structural relaxation time (r)oi salol from photon correlation measurements at different temperatures M (T=267.0 K ([>), 
268.5 K (o), 271.1 K (a), 274.5 K (<)), 278.2 K (©), 280.4 K (□)). d) Structural relaxation time of glycerol from dielectric 
measurements at different temperatures (data taken from Ref. T=217.5 K (<), T=230.6 K (0), T=240.9 K (□), T=247.3 
K (o). In all panels, the solid lines represent the best fit with the PEAG equation Eq. (Q). The relaxation time for dielectric 
measurements is r = l/2iri> rn ax, with Umax the frequency of maximum dielectric loss. The average relaxation time for photon 
correlation data is given by (r) = {tk / /3k )r(/3^ 1 ), where tk and /3k are the Kholrausch- Williams- Watts characteristic 
time and stretching parameter, respectively, and P denotes the Euler gamma function. 



excess entropies, we shall assume that the isobaric con- 
tribution to S c can be expressed, to first approximation, 
as a fraction of the total excess entropy: 

S c (T,0) = *S£?(T,0). (3) 

Hence, we modify the PEAG formula for the structural 
relaxation time as follows: 

C 



t(T, P) = t exp 



T($5, 



ob . 

exc 



Sl soth ) 



(4) 



here, the parameter $ (<1) quantifies the fraction of ex- 
cess entropy at atmospheric pressure configurational in 
nature. Note that the presence of the term S l c soth in 
Eq. (0) prevents the parameter $ to play the role of a 
simple renormalization constant. 

The liquid over crystal excess entropy at atmospheric 
pressure, S\ S °^{T, 0), can be evaluated from calorimet- 
ric measurements. In calculating the isothermal part of 



the configurational entropy, S^, s (T, P), we do not make 
any modification of the expression given in Ref. [^| , since 
no mention to a reference state — crystal or glass — was 
necessary. Accordingly, we approximate the vibrational 
thermal expansivity at P by its value at ambient pres- 
sure, i.e., {dV/8T)f r « (dV/dT) v tbr [Eq. (§)]. This is 
due to the fact that the value and the pressure depen- 
dence of the thermal expansivity for crystals and glasses 
is usually much smaller than that for liquids. Therefore, 
we are left with calculating the integral over pressure of 
the thermal expansivity of the melt. Starting from the 
Tait equation of state ||, simple algebra allows us to 
express this contribution in terms of thermodynamic pa- 
rameters obtained from dilatometric measurements 

Testing the model. — A complete test of the model 
requires the use of information determined by sev- 
eral experimental techniques, i.e, T and P-dependence 
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structural relaxation times, extended calorimetric data, 
pressure-volume-temperature (PVT) measurements. In 
the following we consider three model glass-formers: or- 
thoterphenyl (OTP), phenyl salycilate (salol), and glyc- 
erol. For these systems all the needed data are available, 
and can be properly used to test the generality of our ap- 
proach. Altogether, calorimetric, volumetric and r(T, 0) 
data provide all the parameters in Eq. (EJ) — for details 

f gy \ vibr 



see Ref. 



$ and (§£) 



which can 



— except for C, ^ ...... , ,, , n 

be obtained by a multi-variable fit on experimental re- 
laxation data as a function of T and P. 

Orthoterphenyl. — Isobaric heat capacity of crystalline, 
supercooled and stable liquid OTP determined via 
calorimetry |J has been used to calculate the liquid over 
crystal excess entropy S exc as a function of temperature. 
Specific volume data of the liquid and glassy OTP as a 
function of both temperature and pressure |ll|| have been 
used to determine the parameters of the Tait equation of 
state. Structural relaxation times obtained by means of 
two different techniques are available from the literature. 
A study of pressure effects on the dielectric relaxation of 
supercooled OTP [jl2| reports several isobaric measure- 
ments as a function of temperature which are shown in 
Fig , EJa together with dielectric data at ambient pressure 
|L3| . A second pressure and temperature dependent pho- 
ton correlation investigation has been reported in p[ . 
The average relaxation time (r) appears in Fig. [j]b in the 
T and P ranges investigated. 

We fit according to Eq. (EJ) the data in the T range 247- 
272 K at atmospheric pressure, where log 10 t is linear vs. 
(TSexc)^ 1 ), and all the isobaric dielectric data at higher 
pressures. The simultaneous fit gives: $ = 0.69 ± 0.05, 



(dV/dT) v Q lbr 



10" 



mol" 1 K" 



C = 



(4.3 ± 0.7) x 

(1.49 ± 0.12) x 10 5 J mol" 1 . It is evident from Fig. |a 
the high quality of the agreement between experimental 
relaxation time data for OTP (symbols) and calculated 
curves (solid lines) using the PEAG model in the form 
of Eq. (Q). As a cross-check of the reliability of the fit 
performed on the dielectric relaxation data, we plot the 
same curves as in Fig. EJb on the PCS data. We find that 
the PCS data are perfectly described by the same fitting 
parameters as the dielectric data, except for a vertical 
shift corresponding to a different value of To connected 
with the specific technique used. It is worth stressing 
that the best-fit is set by a physical value of the vibra- 
tional thermal expansivity: the obtained value compares 
well with ~ 5 x 10 -8 m 3 mol^ 1 K . that is expected for 
an OTP-based solid form of matter [jllj . Reasonably, the 
value of the vibrational thermal expansivity corresponds 
to about 30% of the melt thermal expansivity. We re- 
mark that using the PEAG form with a preset $ = 1 
we find an unphysical value for the vibrational thermal 
expansivity (2.3 x 10~ 9 m 3 mol -1 K" 1 ), about two order 
of magnitude lower than the expansivity of the melt. 

Salol. — Precise determination of the heat capacity of 
salol under atmospheric pressure is due to Oguni |L5] . 
Dilatometric measurements V(T, P) of crystalline and 
liquid salol have recently been carried out [fL6[ permit- 



50 



40 



30 



20 



c 

CD 
C 

o 

Q. 

E 
o 
o 

a. 10 
o 

c 

LU 



" OTP 


I 1 1 








(vibr) . 

exc x t T 


_ o o oo9-° t lJ: 


tT - 

## 

s 

exc 


S -0.69-S 

c exc 




i 





210 



240 



270 300 
T[K] 



330 



FIG. 2: Entropy components in OTP at atmospheric pressure. 
Circles represent the excess entropy of the liquid over the 
crystal, S exc = S melt - S crystal , calculated from experimental 
calorimetric data. The best fit with S exc (T) = Soo — k/T for 
T>247 K, and its extrapolation for T<247 K are shown with 
solid and dash line, respectively. The thick dash-dot line rep- 
resents the configurational entropy evaluated as a fraction of 
the excess entropy, S c = &S e xc- Tft=203.8 K is the tempera- 
ture at which the entropy of the liquid equals the entropy of 
the crystal. 



ting to extract the parameters of the Tait equation of 
state. Isothermal relaxation times measured by photon- 
correlation spectroscopy are available (?]]. The pressure 
variation of (r) is shown in Fig. [l]c in the temperature 
range investigated. We analyze the relaxation times data 
using Eq. (§, and find $ = 0.68 ± 0.08, (dV/dT) v Q lbr = 
(3.8±0.7)xl0" 8 m 3 mol" 1 K _1 , and C = (1.9±0.3) x 10 5 
J mol -1 . The value for the vibrational thermal expan- 
sivity is in remarkable agreement with that found for the 
poly-crystal via PVT measurements [{dV / dT) C Q VStal ~ 
(4.0 ±0.5) x 10~ 8 m 3 mol" 1 K" 1 ] gfjj. A previous anal- 
ysis of these data performed by some of us Q| using 
a preset $ = 1, i.e. replacing S c with S exci provided 
{dV/dT)l lbr = (1.0 ± 0.7) x 10- 8 m 3 mok 1 K" 1 , which 
is in feeble agreement with the experimental value of the 
poly-crystal. 

Glycerol. — Isobaric heat capacity of liquid-glass and 
crystal glycerol JlTj have been used to evaluate the liquid 
over crystal excess entrop y. V olumetric measurements of 
glycerol reported in Ref. jL8| as a function of both T and 
P have been represented in terms of the Tait equation. 
As far as r(T, P) data, we focus on the pioneering work 
by Johari and Whalley p9| , which provides relaxation 
time data in agreement with those recently reported by 
Lunkenheimer at atmospheric pressure plj. For our test, 
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we subject to the fit procedure data up to ~ 1 GPa (see 
Fig. [j]d) excluding all isotherms at temperature higher 
than 247 K, in order to remain in the validity range of 
the PEAG model. Because of the not very high qual- 
ity of the data, and with the background of the analysis 
of OTP and salol in place, we undertake the fit proce- 
dure with [dV/dT) V Q br constrained to a physical value. 
Dilatomctric measurements of glycerol provide an aver- 
age value (1.68 ± 0.01) x 10~ 8 for the thermal expan- 
sivity of the glassy state in the temperature range in- 
vestigated. As a result of a global fit procedure we find 
C = (12.8 ±0.4) x 10 5 J mol" 1 , and $ = 0.71 ±0.01 ||, 
this value being surprisingly close to that found in OTP 
and salol. 

In summary, the fundamental information we gain 
from the above analysis is that at atmospheric pressure, 
about 70% of the excess entropy is of configurational ori- 
gin for each one of the molecular systems here considered. 
These results astonishingly agree with the value found for 
S c /S exc (0.69) in selenium ||, where the configurational 
entropy was calculated from neutron scattering data. In 
addition^ we observe that simulation studies on SPC-E 
water fl2J 
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25\ have shown S c and S exc to be propor- 
tional, with a constant ratio of about 0.77 M. 

In conclusion, we have introduced a general method to 



calculate the configurational contribution to the excess 
entropy for molecular liquids. The general validity of the 
method and the reasonableness of the hypothesis made 
are tested on three model glass formers. The results ob- 
tained include the following points: i) Experiments per- 
formed as a function of T and P are actually able to iden- 
tify the configurational fraction of excess entropy, giving 
a positive answer to our initial question; ii) The value of 
the configurational fraction of excess entropy reconciles 
the results obtained by previous laboratory and compu- 
tational tests; Hi) Configurational entropy is confirmed 
to be a key concept in controlling the slow dynamics of 
glass forming systems. 

This last evidence is also confirmed by the energy 
landscape approach to the glass transition (26), as high- 
lighted by recent simulation studies (see, among others, 

Refs. H El H H S EH)" We believe that more ef- 
forts should be devoted in order to clarify the connection 

among the potential energy landscape features (see, for 

instance, the equation of state proposed in Refs. |3l| , |32| ) 

and experimental results in T and P domains. 

The authors would like to particularly thank Prof. M. 
Oguni for shearing calorimetric data on salol prior to 
publication. 
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